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Abstract. The conventional way to search for long-lived CHArged Massive Particles (CHAMPs) is to 
identify slow (small /3) tracks using delayed time of flight and high ionization energy loss. But at the 7—14 
TeV center of mass energy of the LHC, a CHAMP may be highly boosted (high /3) and therefore look 
more like a minimum ionizing particle, while for high momentum muons (^500 GeV/c) the radiative effect 
dominates energy deposition. This suggests a new strategy to search for CHAMPs at the LHC. Using 
energy deposition from different detector components, we construct a boosted decision tree discriminant 
to separate high momentum CHAMPs from high momentum muons. This method increases substantially 
the CHAMP discovery potential and it can be used to distinguish possible di-CHAMP or CHAMP-muon 
resonance models from di-muon resonance models. We illustrate the new method using a mCMSB model 
and a recently proposed di-CHAMP model and we give updated CHAMP mass limits for these two models 
using the results from a recent CDF CHAMP search. 



1 Introduction 

Many extensions of the standard model (SM) suggest the 
existence of a CHArged Massive Particle (CHAMP) that 
is long-lived [l]-^. The long lifetime is caused by limited 
phase space due to small mass splittings or by small cou- 
plings. If the CHAMP escapes a particle detector before 
decaying it can be regarded as stable. CHAMP searches 
play an important role in constraining SUSY models (see, 
for example, Ref. [o]). Therefore, there is interest in pur- 
suing these searches using early LHC data. 

A number of searches for weakly interacting CHAMPs 
were performed at LEPH 10-13J, yielding 95% C.L. lower 



mass limits of ~100 GeV / c for a stable supersymmet- 
ric partner of SM leptons [m]. Tevatron searches for CH- 
AMPs have been carried out recently by CDF 15 and 
DO [16^. The conventional method to detect CHAMPs is 
to search for slow moving and/or highly ionizing parti- 
cles [It]. This technique has been used for the LEP and 
Tevatron searches, where the CHAMPs are expected to be 
slow due to the low center of mass energy and the expected 
high CHAMP mass (>100 GeV/c^). 

At the energy of the LHC (v^ =7-14 TeV), however, 
CHAMPs can be highly boosted. Therefore, they behave 
more like minimum ionizing particles (MIPs), and con- 
sequently, the conventional search technique may lead to 
significant signal losses. Current experimental efforts at 
the LHC [18,19 focus on slow moving CHAMPs, separat- 
ing them EmQ muons using their larger ionization energy 



loss and time of flight delays. In this paper, we propose a 
way to search for high speed CHAMPs. We use the dif- 
ference in energy deposition in the sub-detectors, as high 
speed (high momentum, p) CHAMPs behave like MIPs 
while high p muons tend to deposit more energy due to 
radiation effects. The goal is to increase the discovery po- 
tential of the LHC experiments using early data. 



2 Models Under Investigation 

We investigate two models for this study: the minimal 
gauge mediated supersymmetry breaking (mOMSB) model 
1 and a resonant di-CHAMP model |6 . In all simulations, 
we use = 10 TeV for the LHC and = 2 TeV for 
the Tevatron. The mCMSB model we explore corresponds 
to the SPS7 benchmark scenario proposed in \20^. In this 
model, the stau (f) is the next-to-lightest supersymmet- 
ric particle (NLSP). We vary A from 31 TeV to 100 TeV, 
with fixed parameters Nmes = 3, tan/3 = 15, /i > 0, Cgrav 
= 10000 and Mmes/^ = 2. The large value of Cgrav re- 
sults in a long-lived stau, while yl = 31 — 100 TeV gives 
a stau mass of 100 GeV/c^ to 308 GeV/c^ We produce 
the mGMSB particle mass spectrum and the decay table 
with the program ISASUGRA [21 , which is used as input to 
22 with all SUSY processes enabled. Figure [l] 
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PYTHIA 6.4 

shows the distribution for some kinematic quantities for 
three different values of the stau mass within this bench- 
mark scenario, at an LHC energy of 10 TeV. 

Some recent models predict di-CHAMP 6 , or CHAMP- 
muon 7 resonances in which the CHAMPs are relatively 
light compared to the resonances, and thus obtain a large 
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Fig. 1. Distributions for some kinematic quantities for the 
SPS7 benchmark scenario at stau masses of 100, 156, and 308 
GeV/c^ for: (a) stau p, (b) stau 77, (c) stau /3, (d) di-stau in- 
variant mass (assuming track has zero mass). All distributions 
are normalized to unit area. 
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Fig. 2. Distributions for some kinematic quantities for the di- 
CHAMP model with CHAMP mass 121, 182 and 302 GeV/c^ 
for: (a) CHAMP p, (b) CHAMP ry, (c) CHAMP /3, (d) di- 
CHAMP invariant mass (assuming track has zero mass). All 
distributions are normalized to unit area. 



boost. We generate di-CHAMP events using CompHEP 
4.5.1 [23 with a model file provided by Kilic, Okui and 
Sundrum based on Ref. 6 . The CompHEP events are 
hadronized using the program Pythia. We vary the Mp 
parameter in the di-CHAMP model from 1.0 TeV/c^ to 
2.5 TeV/c^ yielding CHAMP (K) masses between 121 
GeV/c^ and 302 GeV/c^ The K to p mass ratio is 0.12 in 
the model used and the factorization scale is set to the in- 
variant mass of the di-CHAMP system. Figure [2] shows the 
distribution for some kinematic quantities for this model 
at 10 TeV. 

Both models predict events that may contain two CH- 
AMPs. The di-CHAMP invariant mass distributions are 
shown in Figs, [l] and [2] Because high P CHAMPs resem- 
ble muons in their experimental signatures (for example, 
they are deeply penetrating and have a short travel time 
to the detectors), these events could appear in searches 
for high mass di-/i resonances (as predicted by Z' mod- 
els [24|j27] ) or in the high mass tails of di-/i distributions 
(as in contact interaction models [28l[29 |). 

We define /3 > 0.95 as the high momentum (muon-like) 
region and 0.6 < /3 < 0.8 as the traditional CHAMP re- 
gion. The 0.6-0.8 f3 range is the approximate range of the 
planned LHC experimental searches. Both ATLAS and 
CMS have acceptance for P < 0.6, but the efficiency falls 
rapidly 18,19 below this value. Figure [3] shows the frac- 
tion of events in the two regions for the mGMSB and 
di-CHAMP models at 10 TeV. For both models, more 
events are produced in the high momentum region than in 
the traditional search region for CHAMPs of lower mass. 
Given the expected higher trigger and reconstruction ef- 
ficiency for the high momentum region than for the tra- 
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Fig. 3. The fraction of tracks that fall within /3 > 0.95 and 
0.6 < P < 0.8 for the mGMSB stau model (red) and the di- 
CHAMP model (blue) at 10 TeV LHC energy 



ditional CHAMP search region, we would expect to ob- 
serve more high momentum than low momentum CHA- 
MPs when their masses are low. 

Previous experimental results from LEP and the Teva- 
tron limit the interesting phase space for high momentum 
CHAMPs at the LHC. CDF has set a limit of 10 fb on 
the cross section to produce a single, weakly interacting 
CHAMP. To apply this limit to our benchmark models, 
we must account for the fact that some events may have 
two CHAMP candidates. Therefore, we calculate the ac- 
ceptance, a, for an event to pass the CDF slowly mov- 
ing minimum ionizing particles trigger criteria {pt > 40 
GeV/c, \r]\ < 0.7, 0.4 < /3 < 0.9, and isolation) given one 
or two chances: 



a = 



Ntot 



(1) 
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Table 1. Cascade and pair production cross sections at Teva- 
tron and LHC energies for the SPS7 benchmark scenario in 
mCMSB. The last column shows the acceptance for an event 
to have at least one CHAMP within the CDF analysis selection 
for cascade production at the Tevatron energy. 



Table 2. Production cross sections at Tevatron and LHC en- 
ergies for the di-CHAMP model. The last column shows the 
acceptance for an event to have at least one CHAMP within 
the CDF analysis selection at 2 TeV Tevatron energy. 



where Nfot is the total number of events generated, A^i 
(A^2) is the number of events with one (two) CHAMPs 
within the trigger criteria, ei is the efficiency given by 
Ref. [15] to find a CHAMP given one chance, and €2 is the 
efficiency to find a CHAMP given two chances 30 . 

Figure |4] shows the cross section times acceptance for 
these models as a function of the CHAMP mass. We in- 
terpret this plot as excluding stau masses below ~140 
GeV/c^ and di-CHAMP resonances with CHAMP masses 
below -180 GeV/c^ Tables [l] and [2] list the cross sections 
for the two models, and the corresponding CDF accep- 
tance at ^/s = 2 TeV. Based on these exclusion limits, 
we select the stau mass = 156 GeV/c^ in the mGMSB 
model and the K mass = 182 GeV/c^ in the di-CHAMP 
model as our benchmark signals (labeled staul56 and di- 
CHAMP 182 in the following). 



3 Event Generation and Selection 

In addition to the signal (see Sec. [2|, we model multiple 
sources of SM (muon) backgrounds (W, ^/7*, multijet, 
tf, and diboson) using Pythia 6.4. All samples are gen- 
erated using the CTEQ 5L parton distribution functions 
(PDFs) [31]. We will use CMS as our reference LHC de- 
tector for this study. However, the principles (modified 
appropriately) are also applicable to ATLAS. 

To select CHAMP candidate events, we require at least 
one high momentum {p > 500 GeV/c) muon-like track 
within the |?7| < 1.479 region. The 500 GeV/c momentum 
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Fig. 4. Cross section times acceptance for the mGMSB SPS7 
benchmark scenario (red) and di-CHAMP model (blue) at 
y/s — 2 TeV. The horizontal black line shows the CDF limit 
on single, weakly-interacting CHAMPs ^15j . 
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requirement corresponds to (3 > 0.95 for the staul56 and 
-(3 > 0.94 for di-CHAMP182, which are fast enough to ap- 
pear similar to muons (for masses at the upper edge of the 
range considered in these models, this requirement corre- 
sponds to /3 > 0.85). Moreover, this ensures that back- 
ground muons are highly boosted and their energy loss 
is dominated by radiation. The muon critical energy oc- 
curs at a muon momentum of 317 GeV/c in copper, 170 
GeV/c in lead tungstate and 581 GeV/c in silicon |32 . The 
|7^| < 1.479 selection corresponds to the barrel part of the 
CMS electromagnetic calorimeter (Ecal) |33 , where the 
energy measurement is most accurate. The rj requirement 
also helps reject high momentum muon background in the 
forward region. In addition, we require a second muon- 
like track with momentum p > 100 GeV/c and pt > 20 
GeV/c in the I77I < 2.5 region to reduce the iy+(jets) 
background. Table |3] shows the cross sections for the sig- 
nal and background processes after the leading and second 
leading CHAMP/muon requirements. 

After these simple selection criteria, the background 
level is already much smaller than, or similar to, the ex- 
pected signal for the staul56 and di-CHAMP182 bench- 
mark models. Therefore, even a simple counting experi- 
ment may be sensitive to a CHAMP signal. However, due 
to PDF uncertainties and the fact that Pythia is only 
a tree level generator, the background assumed here may 
be underestimated. The effect of a larger background will 
be discussed later. In any case, the energy deposition in 
the detector can be used to separate further signal from 
background. 



4 Detector Simulation 

We use Geant4 |34| to simulate a simplified model of 
the CMS detector f33l, shown in Fig.jsJ In our model, the 
tracker is made of 20 layers of 300 jam silicon, the electro- 
magnetic calorimeter is a 5 by 5 array of PbW04 crystal 
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Table 3. Cross sections for signal and background samples be- 
fore selection and after leading and second leading candidate 
selections. The VK/Z+(jets) and multijet backgrounds have 
generator level cuts of 30 GeV/c and 215 GeV/c on the leading 
jet pt, respectively. 
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Fig. 6. Distributions of (a) dE/dx in the Tracker, (b) energy 
deposition in Ecal, and (c) energy deposition in Heal. The high 
momentum staul56 signal and the background are shown with 
a solid and an empty histogram respectively. 



Fig. 5. GEANT4 simulation: a muon track with p = 650 GeV/c 
passing through a simplified model of the CMS detector. 



with the same dimensions as a CMS Ecal crystal, and 
the hadron calorimeter (Heal) is composed by 30 layers of 
brass with transverse size similar to that of a CMS Heal 
tower. No magnetic field is included in our simulation, as 
we expect high momentum tracks to be rather straight in 
the detector. 

The energy loss per path length (dE/dx) for a muon or 
a CHAMP in the tracker is calculated using the truncated 
mean f35^ estimator in which the arithmetic mean is com- 
puted after truncation of the highest 40% of the charge 
samples for each track. We take the Ecal energy to be the 
energy deposited in the crystal traversed by the CHAMP 
or muon. This was chosen to minimize the possible en- 
ergy contribution from nearby charged or neutral particles 
and to avoid the need to simulate the shower shape in a 
magnetic field. The Heal energy is computed by summing 
the energy measurements for all layers of Heal, which is 
roughly the same as the single tower energy defined in the 
CMS detector. 

To better account for the detector response, the Geant4 
simulated energy is smeared with the CMS detector energy 
resolution functions 33 . For Heal it is cr/£: = 120%/ 
6.9% 36 . The Ecal energy is smeared using {a jE)^ — 
(2.8%/V^)2 + (0.12/£^)2 -h (0.30%)2. For tracker dE/dx, 
we do not smear the energy because the simulated energy 
resolution is approximately the same as that of the re- 
constructed energy. This is due to the very high signal to 
noise ratio for the CMS tracker readout electronics [33^. 



As noted in Secjsj the main background is Z/7*+(jets) 
after our simple selections. We save all high momentum 
tracks from all background or signal events and create 
background samples by mixing events according to their 
cross sections. The muon/CHAMP tracks are then passed 
to our simplified CMS GEANt4 simulation and we obtain 
the energy deposition for each track. 

The results of our simulation and energy smearing in- 
dicate that the distributions of CHAMP energies in Ecal 
and Heal should be Gaussian (see Fig. [6|. This is due 
to the fact that CHAMPs primarily lose energy through 
ionization. For the thick material in Ecal and Heal, the 
Landau distribution becomes Gaussian. The high speed 
muons, however, exhibit long tails in dE/dx distributions, 
corresponding to radiative energy losses. For the tracker, 
the dE/dx distribution is Gaussian in both cases, as ex- 
pected for the truncated mean estimator 35 . 

In Fig.[7| we show the average values of Tracker dE/dx, 
energy deposition in Ecal, energy deposition in Heal, and 
energy deposition difference between Heal and Ecal as a 
function of momentum, for two CHAMP models and two 
mass points. The low momentum region (< 500 GeV/c, 
corresponding to low /3s, which is the focus of conventional 
CHAMP searches) shows the high ionization energy loss, 
which is inversely proportional to ff^ |32 . We make an im- 
portant observation: the average energy deposition in Ecal 
and Heal in the high momentum region (> 500 GeV/c) are 
almost identical for the four model/mass benchmarks con- 
sidered. Note that the muon energy deposition increases 
with momentum, as expected from the increased radiative 
contributions, while the Tracker dE/dx of muons is evenly 
distributed in approximately the full momentum range, 
mainly due to the fact that here dE/dx is calculated with 
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Fig. 7. The average (a) Tracker dE/dx, (b) energy depo- 
sition in Ecal, (c) energy deposition in Heal, and (d) differ- 
ence between Heal and Ecal in energy deposition, as a func- 
tion of CHAMP/muon momentum for four different CHAMP 
model/mass points and background muons. Further analysis 
in this paper will concentrate on the high momentum (> 500 
GeV/c) region. 



40% truncation, which does not take into account high 
energy (radiative) hits. 

Given the important observation that the energy de- 
positions in these detectors is largely model and mass in- 
dependent in the high p region (Fig. [t]), we can use the 
energies to develop a new method to separate CHAMPs 
from muons in a manner that is analogous to standard 
particle identification that separates electrons and pho- 
tons from hadrons and jets |37 . In the following section 
we use a multivariate method to enhance the discrimina- 
tion between muons and CHAMPs. 




0.4 0.6 0.8 
BDT Discrimirtant 

Fig. 8. Distributions of the BDT discriminant for staul56 
signal (solid) and background (empty). 



5 Multivariate Analysis 

The direct use of the energy depositions, whose distribu- 
tions are shown in Fig. [6) can only provide limited sep- 
aration between CHAMPs and muons passing our basic 
selection. However, it is well established that better sep- 
aration between distributions can often be obtained by 
combining well understood variables into a multivariate 
discriminant. To that end, we use the boosted decision 
tree (BDT) method provided in the Toolkit for Multivari- 
ate Data Analysis (TMVA) 38 to compute discriminants 
that take the three energy depositions as inputs. 

Figure [8] shows the BDT discriminant distribution for 
the staul56 model and backgrounds after our simple selec- 
tion. We see very clear separation between high speed CH- 
AMPs and background muons. The di-CHAMP182 model 
has a behavior very similar to the one of staul56 model. 
We estimate the significance of the signal using 



S,L = \l2{{S + B)\n{l + ^)-S) 



(2) 



where S and B are the number of signal and background 
events fSO] that remain after a cut on the BDT discrimi- 
nant. The cut is chosen by maximizing Scl- 

We use ScL to illustrate the expected improvement in 
sensitivity of the BDT method relative to the conventional 
one. A more precise estimate of the sensitivity of the LHC 
experiments for observing a possible CHAMP signal would 
require a full detector simulation. However, the main point 
we wish to make is that a better use of the energy deposi- 
tions can enhance significantly the discovery potential of 
the LHC using the early data. The quantity ScL is suffi- 
cient to make this point. 



6 Application of the BDT Discriminant 

6.1 Discovery potential for CHAMPs 

For the two benchmark models. Table [4] compares the 
ScL values expected using a simple counting method with 
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Table 4. Signal significance assuming 200 pb ^ data at 
10 TeV LHC, using simple counting method or BDT. 





nominal 


double nominal 


model 


background 


background 




counting 


BDT 


counting 


BDT 


staul56 


5.6cr 


8.3cr 


4.4cr 


7.1a 


di-CHAMP182 


4.2cr 


6.7cr 


3.2cr 


5.8cr 



those expected from the BDT method. We have assumed 
an integrated luminosity of 200 pb~^, collected at an LHC 
center of mass energy of 10 TeV. We note that a simple 
counting experiment can already give an ScL as high as 
5.6(7 (staul56 model with nominal background) using our 
selection of high momentum tracks within the muon ac- 
ceptance. Even without looking at the energy deposition 
differences between muons and CHAMPs, a low mass CH- 
AMP discovery is already possible for some model/mass 
points. 

The background could be larger than our estimate due 
to several effects. The PDFs are not well constrained at 
LHC energies. Our simplified GEANT4 simulation does 
not include track momentum mis- measurement that could 
cause moderate momentum (100-500 GeV) muons to re- 
construct at larger momentum (> 500 GeV). A full de- 
tector simulation and reconstruction as well as compari- 
son to real data is necessary to incorporate these effects. 
However, Table [4] also shows the effect of doubling the 
background estimate to simulate these types of effects. 
The larger background would lower the significance of 
both methods. However, the BDT discriminant can en- 
sure discoveries (> 5cr significance) for the staul56 and 
di-CHAMP182 benchmarks, even if the background were 
larger than our nominal estimate by a factor of two. 

The method proposed here has very high detection effi- 
ciency because the CHAMP is moving very fast and there- 
fore, like a high pt muon, its trigger and reconstruction 
efficiencies are very high at CMS and ATLAS. (For P less 
than ~0.6, however, the trigger and reconstruction effi- 
ciencies for CHAMPs drop significantly at CMS 18 and 
ATLAS [19].) Moreover, the new method can use high mo- 
mentum cosmic muons to evaluate instrumental effects re- 
lated to energy depositions. Recent results from the CMS 
collaboration show the GEANT4 simulation is consistent 



with detector response to cosmic muons 40 -42 . 

We then compare the expected significance of the low 
/3 method, the high p counting method, and the new BDT 
method for the two benchmark models. Table [5] lists the 
ScL values from the three methods, assuming 200 pb~^ 
data from 10 TeV running at the LHC, for several mass 
points. The efficiency is obtained from the number of CH- 
AMP tracks in the slow (0.6 < (3 < 0.8 and \r]\ < 0.8 and 
Pt > 40 GeV/c) or fast CHAMP regions (as defined in this 
paper), divided by the total number of events. We assume 
100% trigger and reconstruction efficiency for CHAMPs in 
both the conventional and BDT methods, a background 
of < 1 event for the conventional method 18 , and our 



Table 5. Signal significance from the conventional, counting, 
and BDT method, assuming 200 pb~^ data from 10 TeV LHC. 
The efficiency for the events to fall within the basic selection 
(see text) is given in parenthesis. 
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Fig. 9. Significance of the new and conventional search meth- 
ods as a function of the CHAMP mass for the two models 
explored: mGMSB SPS7 (red) and di-CHAMP (blue). 



Figure [9] shows the significance of the new and con- 
ventional methods as a function of CHAMP mass. We see 
that for the masses of interest for early LHC running, the 
high momentum analysis has comparable or better dis- 
covery potential. The curve for the conventional method 
is a rough estimate and it could turn out to be larger, 
but we do not expect it to be radically different. The two 
search methods are complementary as they explore differ- 
ent regions of signal acceptance. Therefore, the discovery 
potential could be enhanced further by combining the two 
methods, taking due account of correlations. 

Although this paper uses CMS as a model detector, 
the method is applicable to ATLAS, which has similar 
Heal and Ecal energy resolutions as CMS [43]. Although 
ATLAS does not measure the Tracker dE/dx in the same 
way as CMS, the number of hits in the Transition Radia- 
tion Tracker (TRT) will give similar discriminating power. 
This is because a high speed muon results in more high 
threshold hits, but a high speed CHAMP (still in the (3^ = 
1 — 10 region) only gives about 2 hits, as pointed out in 
Ref. 144. In addition, the technique could be extended to 
include information from other ATLAS or CMS detectors 
such as proposed in Ref. [45]. 
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Fig. 10. Di-/i mass for main background and signals (staul56 
and di-CHAMP182 models, unstacked) after initial selections 
in this analysis. The di-/i Mass is smeared with mass resolution 
in Ref. 46 . Plot is scaled to 200 pb'^ 10 TeV LHC data. 



goal. By combining the two search methods, one can, in 
principle, increase the discovery potential further. Finally, 
the BDT discriminant provides an experimental tool to 
distinguish CHAMP models (di-CHAMP, CHAMP-muon 
resonance models or non-resonance CHAMP models) from 
models with di-muon excesses in standard high mass di-/i 
searches. 
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6.2 Distinguishing di-/i signal from di-CHAMP or 
CHAMP-/i signal 

In the di-CHAMP model, CHAMPs can originate from the 
decay of a heavy resonance. Therefore, it is interesting to 
look at the invariant mass spectrum (calculated assuming 
the muon mass) for those events. As shown in Fig. 10 , after 
the initial selection used in this analysis (not including the 
BDT), we see a broad excess above 500 GeV/c^ in the di-/i 
mass spectrum. Previous Monte Carlo studies have shown 
some evidence that a peak may be observable [47] , but it 
is not seen here with this benchmark model. 

Even for CHAMPs that do not originate from a com- 
mon resonance, both models predict an excess of events 
at high di-/i mass (Fig.[lO|. Therefore, an analysis search- 
ing the di-fi spectrum above the Z peak (for evidence 
of a Z', compositeness, Randall-Sundrum gravitons, etc.) 
will need to determine whether any observed excess arises 
from muons or from CHAMPs. Our new method can eas- 
ily distinguish events with two CHAMPs from those with 
two muons. The method will also be useful in identify- 
ing CHAMP-muon resonances such as those suggested in 
Ref. [7]. 



7 Conclusion 

We propose a way to search for high speed CHAMPs at the 
LHC experiments to complement the conventional slowly 
moving CHAMP searches. The SM background can be 
kept to a very low level using simple selection criteria. 
With 200 pb-^ data collected at 10 TeV during the first 
run of the LHC, a simple counting experiment already 
has the possibility to discover a stau with mass just above 
140 GeV/c^. To improve the signal-background separa- 
tion, we construct a boosted decision tree whose inputs 
are the energy depositions from muon-like tracks. This 
discriminant can improve significantly the discovery po- 
tential during early LHC running, which is our principal 
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